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1. INTRODUCTION 


Dr. Gautam Badhwar of the National Aeronautics and Space Administration (NASA), 
Lyndon B. Johnson Space Center (JSC) has developed an approach to classifi- 
cation based on the hypothesis that a given crop, after emergence, has a 
unique spectral profile in time in the spectral region covered by the Land 
Satellite (Landsat) channels. A crop exhibits emergence at different times - 
over a site. The result is displacement of the crop profile, but the profile 
retains the same shape (ref. 1), This classification method incorporates the 
effects of emergence data distribution and bases classification on the tem- 
poral profile of the crop (ref. 2). Local training, one input training field, 

l 

is us^d to establish the crop profile for the site. 

The classification maps of 26 corn/soybean sites were produced using this 
method for evaluation as an aid to quality assurance for the Accuracy Assess- 
ment {AA) digitized ground-truth inventory maps (ref. 3). Since this method 
is in the development stage, it was expected that the generation of results 
for so large a data set would (]) test the applicability of the program to 
varying acquisition distributions and training field signatures and (2) define 
a tentative method for producing an acceptable classification as well as an 
initial set of criteria for evaluating the acceptability of the classification 
maps. 

In this report, the software and the procedures used to apply Badhwar profile 
comparison classification are documented, and numerical results are presented. 
Section 2 defines the data set: site selection, acquisition choice, and 
editing of picture elements (pixels) using "SCREEN." Section 3 describes the 
available software programs. Section 4 discusses two aspects of the analysis: 
the analysis necessary to define the input to the program and the analysis 
u$ed to evaluate the classification output. Section 5 presents numerical 
results in tabular form. Conclusions and recommendations, including a recom- 
mended data set for further program development, are given in section 6. 
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Mr. Gary L. Gutchewski of NASA/JSC selected the training fields for this study. 



2. DATA SET 


Twenty-six corn/soybean segments were classified by using acquisitions for 
which the distribution over the corn-growing season was adequate to define a 
curve (crop profile) in each Landsat channel. The full segment (22 932 pix- 
els), less those pixels removed by the SCREEN program developed by the 
Environmental Research Institute of Michigan (ERIM), was classified. 

2.1 SITE SELECTION 

Twenty-nine segments in the U.S. Corn Belt were used for the Multicrop 
Exploratory Experiment Test (ref. 4). These sites were located in agrophysi- 
cal unit (APU) 14, 24, 25, or 28, had reasonable acquisition histories, and 
had available the aircraft photography ground truth. Final site selection 
for the multicrop experiment was determined by the procedural constraints on 
the acquisition requirements; of the 29 segments, procedure 1 (PI) estimates 
were made on 25 segments. Segments excluded from the test were sample seg- 
ments 216, 878, 891, and 893 [sample segments 202 and 807 were not used in 
the followup Simulated Aggregation Test (ref. 5) because of poor data quality 
in crucial acquisitions]. 


For use of Badhwar classification, acquisition requirements are different 
from those for analyst labeling and for PI classification. Hence, the entire 
set of 29 segments proposed for the multi crop test .was taken as the basic 
site data set: 

Sample • * . Sample 


segment 

County, state 

APU 

segment 

County, state 

APU 

107 

Boone, 111. 

25 

202 

Atchison, Mo. 

14 

123 

Hamilton, Ind. 

28 

205 

Clark, Mo. 

25 

127 

Montgomery, Ind. 

28 

216 

• Mercer, Mo. 

25 

133 

Whitley, Ind. 

28 

800 

Clinton, Iowa 

25 

135 

Chickasaw, Iowa 

24 

807 

Henry, 111. 

25 

141 

Madison, Iowa 

25 

809 

Ogle, 111. 

25 

144 

Wapello, Iowa 

25 

832 

Adams, Ind. 

28 
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Sample Sample 


segment 

County, state 

APU 

segment 

County, state 

APU 

837 

Benton, Ind. 

28 

877 

Ida, Iowa 

14 

842 

Henry, Ind. 

28 

878 

Kossuth, Iowa 

24 

843 

Henry, Ind. 

28 

880 

Monona, Iowa 

14 

852 

Randolph, Ind. 

28 

881 

Monona, Iowa 

14 

853 

Randolph, Inc. 

28 

882 

Palo Alto, Iowa 

24 

860 

Wells, Ind. 

28 

891 

Shelby, Iowa 

14 

864 

Crawford, Iowa 

14 

893 

Webster, Iowa 

24 

865 

Crawford, Iowa 

14 





2.2 ACQUISITION SELECTION 

The Badhwar classification method reported in this document requires five acqui- 
sitions in the .postemergence to preharvest growth stages of corn. Classifica- 
tion is successful, however, if four acquisitions are available in this period. 
Site selection depends upon the availability of a suitable set of acquisitions. 

e Acquisition distribution must be adequate to define the crop profile. 

For sample segment 107, available acquisitions 

sampled only the ripe to senescent crop stages 

(Julian days 208 to 305); hence, a proper crop 
profile could not be defined, and the segment 
could not be classified. 

• Acquisitions on cloudy and hazy days should be avoided. The Badhwar pro- 
gram is not, however, overly sensitive to atmospheric conditions; cloud 
shadows and .haze apparent on the film products of the Large Area Crop 
Inventory Experiment (LACIE) may not affect profiles or classifications. 

Acquisitions available for sample segment 807 
are theseV 78137 (popcorn clouds), 78164 (haze 
and data drop), 78209 (haze obscured), 78218 
(one-third haze obscured), 78272 (harvest 
beginning, on corn), 78181, and 78290. This seg- 
ment could not be classified because of poor 
data quality in the preharvest acquisitions. 

• Acquisitions should occur in the growing season of the crop, although con- 
siderable robustness was exhibited by the present technique in the 



definitions of postemergence and preharvest. As with atmospheric condi- 
tion effects, final judgment of the usable acquisitions was made with 
reference to graphs of the proposed training fields over the available 
acquisitions. 

The acquisitions available for sample segment 
893 (78131, 78221, 78266, and 78293) clearly 
did not satisfy the constraint of a minimum of 
four acquisitions in the postemergence to pre- 
harvest period. This segment could not be 
classified. 

2.3 USE OF SCREEN 

As each image was unloaded from an Earth Resources Interactive Processing 

System (ERIPS) image unload tape onto a disk for processing on the programmed 

data processor (PDP 11/45), it was edited using the ERIM program SCREEN 

(ref. 6) — a "procedure for automatically detecting garbled data, clouds, 

2 

snow, cloud shadows, and water in MSS data." Pixels failing to pass this 
edit step were excluded from processing. Since the number of pixels removed 
by SCREEN varied with segment and acquisition, percentages were normalized 
to 22 932 pixels for all segments for the presentation of results in section 5. 


2 . > 

A listing that tallied the number of pixels edited from each acquisition was 
available. This report affected acquisition choice in that overly screened 
acquisitions were avoided, if possible. If it was not possible to avoid 
using an acquisition. with an excessive number of pixels removed, a program 
was available which could be used to omit an acquisition from editing 
(refs. 7 and 8). This was done on all acquisitions for sample segment 135 
(agricultural land was removed by SCREEN as being cloud shadows) and on 
acquisition 78165 for sample segment 144 (10 858 pixels were removed by 
SCREEN as clouds although the acquisition is cloud free). 
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3. SOFTWARE AVAILABLE 


All data processing used to generate the classifications was done on the 
PDP-11/45 image processor. Several software programs were used: IMUNLD, 

IMAPLT (to help define the analyst input training field and acquisition set), 
CLASFY, A2SGMAP, and MISMAP (to aid in analyst evaluation of results). This 
section describes each of these programs. 

3.1 IMUNLD 

IMUNLD takes an image unload tape generated on the Earth Resources Interactive 
Processing System (ERIPS), edits it using SCREEN, adjusts the Landsat-3 acqui- 
sitions into a data range comparable to the data range of Landsat-2 acquisi- 
tions using the Wehmanen multiplicative factors (ref. 9), and loads the images 
into a PDP-11/45 disk. 

• Input: ERIPS image unload tape. 

• Output: screened, Landsat-3 adjusted images on a PDP-11/45 disk. 

3.2 IMAPLT 

3 

IMAPLT (ref. 10) plots the individual pixels of a field, giving reflectance 
values versus time (i.e., the acquisition dates specified) for each channel. 
IMAPLT then plots the field mean values, each channel, with a one standard 
deviation envelope; a curve is fitted through the mean values. Eight graphs 
are produced for a field over a set of acquisitions. Graphs are displayed on 
the Image-100 Tektronix screen, and hardcopies are made automatically. The 
segment number, acquisitions used, coordinates of the field, channel number, 
the number of pixels in the field, and the mean and standard deviation on each 
acquisition are listed on the first plot. The constant values computed from 
the data for the model (with the estimated error), the estimated field plant- 
ing date (with error), the values of the fitted curves at the specified acqui- 
sitions (which can be compared with the computed mean values of the data). 


3 

Available reference is to TRJPLT, an early version of IMAPLT. 


5 



and the chi-square value for the fit of the approximating curve to the field 
data are presented on the second plot. 

• Input: field coordinates (line, pixel) in order; acquisition set of four 

or of five acquisitions. 

• Output: eight graphs as above, two for each Landsat channel. 

3.3 CLASFY 

CLASFY (ref. 11) computes the constants for the curves from the training field 
data, compares (with this crop profile in each channel) the values for each 
pixel in the segment, and rejects those pixels which are not within a speci- 
fied chi-square measure of the profile. The technique for rejection is to 
compare pixel channel values with the profiles in channel 2, channel 3, chan- 
nel 4, then channel 1 in succession and reject if the comparison in any single 
channel is inadequate. Variability of the time of planting/emergence is 
allowed for in the comparison of individual pixels with the crop profile 
(refs. 1 and 2). Accepted pixels are labeled "corn"; rejected pixels, non-corn. 

• Input: five, or four, image files; coordinates of one crop-of-interest 

field to establish crop profiles; initial values for the function constants 
as computed in IMAPLT (to aid convergence of the approximating curve). 

• Output: classification file on disk which has a designation of "screened," 

"corn," or "non-corn" for each pixel in the segment; lineprinter sheet 
summarizing the following: 

a. Acquisitions used. 

b. Training field coordinates and the number of pixels in the field. 

c. Mean and standard deviation for each channel and each acquisition 

i 

(field averages) .. 

d. The input and the final constants (with error) for the model. 

e. Final chi-square values for each channel (training field data). 

f. Estimated planting date of the training field (with error) as derived 
for each channel. 
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g. The chi-square thresholds in each channel applied as cutoff values in 
classification. 

h. The number of pixels cut for exceeding the chi-square threshold, hence 
removed from consideration as corn, in each channel. 

1. The final numerical results: the number of pixels classified as corn, 

the number of pixels screened, and the number of pixels rejected as 
corn. 

3.4 A2SGMAP 

A2SGMAP provides a full classification map (22 932 pixels) of the results of 
CLASFY . The scale is the same as that used for the AA digitized ground- truth 
maps. Pixels classified as corn are designated "C,“ those screened are "T," 
and those rejected as corn are left as blank spaces on the map. 

• Input: classification file from CLASFY. 

» Output: lineprinter map, full classification. 

3.5 MISMAP 

MISMAP (ref. 11) compares the classification map presented by A2SGMAP with the 
AA digitized ground- truth inventory map for the segment. A lineprinter map 
with this code is generated: 

• Pixels screened appear as T. 

• Ground- truth corn classified as corn appears as C. 

• Ground- truth non-corn rejected as corn is left blank. 

0 Ground-truth non-corn classified as corn appears as +. 

0 Ground-truth corn rejected as corn appears as -. 

0 Pixels for which ground truth is not available but which are classified 
as corn appear as $. 

0 Pixels for which ground truth is not available but which are rejected 
as corn appear as %. 

A numerical scene summary is given in confusion matrix form. 
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MISMAP maps can be generated for all pixels or for pure (AA definition) pixel 
only. Pure pixels (AA) are those which on a subpixel level contain only one 
crop. 

• Input: classification file from CLASFY and ground-truth inventory map 

file. 

t Output: full scene lineprinter map comparing the classification map with 

the ground- truth map and a confusion matrix numerical summary of results. 



4. DISCUSSION OF ANALYSIS 


Use of Badhwar classification involves different analyst requirements from 
those necessary to implement a PI classification (ref. 12). Four differences 
are as follows: 

• Badhwar classification is trained on one class — the crop of interest (in 
this case, corn) as opposed to PI, in which training is done on both the 
crop and the non-crop classes. 

• Input for Badhwar classification is one analyst-selected field as opposed 
to PI analyst-labeled dots from a semi predetermined set of pixels that 
samples the major crop categories in the scene. 

9 In Badhwar classification, four or five acquisitions are needed to deter- 
mine the crop profile; these, preferably, are well distributed in the post- 
emergence to preharvest period. In PI , a set of acquisitions, a maximum 
of four, is used; it includes, preferably, a pre-emergence and a postharvest 
acquisition. 

9 Output from Badhwar classification is a full classification map that must 
be evaluated as opposed to that of PI, the monitoring of a pixel-based bias 
correction applied to a statistical estimate. 

Discussion of the analyst input and analysis will be divided into two parts: 

4.1 -- training field and acquisition selection and 4.2 - evaluation of clas- 
sification results. 

4.1 ANALYST INPUT TO CLASSIFICATION: TRAINING FIELD SELECTION AND CHOICE 

OF ACQUISITIW ‘ “ 

Initial assessment of the acquisition coverage for each segment was done by 
the analyst using the LACIE film products. For each segment, acquisitions 
available in the corn growth period were listed and comments presented on data 
quality and crop growth stage for each acquisition. Based on the list and the 
comments, a recommendation of "possible" or "reject" was given for the segment. 
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Criteria for "possible" acquisition coverage of a segment were as follows: 

• A minimum of four acquisitions, reasonably haze-free and cloud-free, in 
the postemergence to preharvest growth stages of corn was necessary. 

• Acquisitions were within LACIE acceptable registration error and were of 
reasonably good data quality. 

For those segments with "possible" acquisition coverage, training fields of 
corn were defined without reference to ground truth. Since fields -of corn 
vary in planting date and development, it was expected that the acquisitions 
available would characterize the crop profile with varying results between 
training fields. Approximately four candidate training fields were selected 
for each segment, and final field choice was determined by the profiles 
generated. 

Criteria for training field selection from imagery are these: 
t Training field size of 20 to 40 pixels. 

• Exclusion of border and edge pixels. 

• Avoidance of roads, drainage patterns, etc., in fields, if possible. 

» No selection of fields with unusual signature (e.g., only irrigated field 
in dryland area). 

• Training field cloud-free and haze-free on all acquisitions used, not 
harvested on final acquisition, and (preferably) some signs of emergence 
on first acquisition. 

The segments with fields defined by the analyst are those in which the pro- 
gram could be expected to be successful. These .segments, indicated by an 
asterisk(*) in table 1, are the following sample segments: 


107 

837 

864 

127 

842 

865 

133 

843 

877 

141 

852 

881 

216 

853 

882 

809 

860 



in 



Processing was attempted for the remaining segments to “test the limits" of 
the program; the author defined two training fields, using ground truth, for 
these segments. Results for the following sample segments should be con- 
sidered exploratory: 


12 3 

807 

135 

832 

144 

878 

202 

880 

205 

891 

800 

893 


All training fields were graphed using IMAPLT. Five acquisitions were graphed 
if available. If a five-acquisition set was not available or if all training 
fields failed criteria using five acquisitions, four acquisitions were graphed. 

Acquisition selection criteria are listed below: 

a. Distribution of acquisitions relative to the training field crop develop- 
ment must be adequate to characterize the crop profile. As mentioned 
previously, it was expected that the training fields would differ in 
planting date and development. Hence, over a given set of acquisitions, 
the crop profile definition would differ between training fields. Fig- 
ure 1 illustrates channel 2, training field 1, sample segment 882 (acqui- 
sitions 157, 186, 222, 231, 267); figure 2 shows training field 2 (same 
segment and same acquisitions). In development, training field 2 appears 
to be slightly later; for these acquisitions, field 1 is approximated 
better by the curve. Both of these fields are acceptable curves. Fig- 
ure 3 illustrates channel 2 "curve" definition for sample segment 107 over 
available acquisitions (208, 226, 235, 244, 262); this definition is un- 
acceptable as the field is "early" relative to all available acquisitions. 
Figures 4 and 5 illustrate an evident difference in crop development. Fig- 
ure 5 shows a field which is later than that in figure 4 in exhibiting the 
characteristic dip associated with chlorophyll absorption. Both figures 
are for sample segment 860, channel 2 profile (acquisitions 160, 197, 232, 
261? 268). Using these acquisitions, the curve generated for training 
field 2 (fig. 5) is not an adequate approximation of the data; training 
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SEGMENT- e88S NO.-PIXS- S3 CHS ACG OATES 159 186 2S3 231 36? 

FIELD CORDS ARE* 35. 93. 34. 98. 38. 99. 39. 93. 



CALENDAR DAY 


PLOT MEANS AND +✓- ONE STAND DEO 
3 4883 >9.8385 -1.8844 1.4578 



CALENDAR DAY 

Figure 1 Sample segment 882, training field 1. 
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SEGMENT-0883 NO -PIXS- 40 CHS ACQ DATES i59 ■ 
FIELD CORDS ARE- 77. 77. 76. 84. 81. 81 

29 3 4.8 12.9 0.9 11.6 0.9 11.1 0.9 16.0 


40.00 


30.00 


20.00 


10.00 


3 6467 -9.6 

25 6752 1,6 

fitted values Ar 


40.00 


30.00 


20.00 


10.00 


Figure 2 ; — . 




SEGMENT-0107 NO.-PIXS- 38 CHS ACQ DATES 203 226 235 244 268 

FIELD CORDS ARE 1 51 128. 51. 135. 58. 133. 58. 131. 



CALENDAR DAY 


PLOT MEANS AND +/- ONE STAND DEU 
3 3964 -0.0378 -0 1171 0.0033 



CALENDAR DAY 


Figure 3.— Sample segment 107, training field 3. 
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seonENT-esee no -pixs- s ? ch2 acq dates 160 is? E32 asi 268 
FIELD CORDS ARE- ??. 144. 77. 149. 81. 1S1 . 82 147. 

27 8 2 5 15,7 1.8 11.4 1.0 13.9 0.9 13.0 1.0 


40 00 

30.00 

20.00 

10 . 00 ' 

0 . 00 ~ 1 1 I ' I 1 i 1 1 r 

.140.00 168.00 196.00 224.00 252.00 

CALENDAR DAY 



3 6609 
26.1821 1 5140 

fitted values Are 


PLOT MEANS AND +✓- ONE STAND DEV 
-6 2233 -0.5415 1.4590 

0.1690 9.7494 

27.6? 15.50 12.63 12.72 13.65 

CHISQR- 0.3855E+01 


40 00 


30.00 


20.00 


10.00 



0.00 


r 


140.00 168.00 196.00 224.00 252.00 


CALENDAR DAY 


Figure 4.— Sample segment 860, training field 1. 
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field 2 is not acceptable as a basis for classification. The data in 
training field 1 (fig. 4) is adequately approximated by the curve; train- 
ing field 1 is acceptable for use in classification. By proper choice of 
a training field and an acquisition set, a representative crop profile was 
defined in each channel for use in classification. 

b. Chi-square fit values in each channel should be less than 10. This cri- 
terion is dependent upon the acquisition selection criteria just des- 
cribed. Notice that the chi-square value in figure 3(b) is 2.5 and 

3.0 in figure 2(b). A misleadingly good chi-square fit is also generated 
when the data are dispersed; many different curves fit well for dispersed 
data with large standard deviations. Chi-square values give an estimation 
of the adequacy of the crop profile curve as an approximation of the train- 
ing field data. It is sufficient only if (1) a curve is defined and (2) the 
data are compact enough that representation by a single curve is sensible. 

c. Estimated planting date for the training field as generated on each chan- 
nel should be the same within the estimated planting date error. This cri- 
terion was seldom violated since the estimated planting date error tended 
to be quite large. 

As mentioned in section 2.2, data quality, atmospheric effects, and crop growth 
stages were assessed by using the film products; but for the Badhwar classifi- 
cation method, final judgment of usability depended upon the curve definition. 
For sample segment 144, day 165, the imagery showed poor data quality. How- 
ever, this is the only date usable for curve definition in early growth stages, 
and it proved to be adequate. Atmospheric effects were also evaluated from 
the graphs. If an apparently hazy, but necessary, acquisition graphed success- 
fully into the crop profile, it was accepted for use. If use of a cloudy 
acquisition could not be avoided, training fields were selected free of clouds; 
and resultant screening-out or misclassification of cloudy corn fields was 
accepted as a penalty. Very early and very late acquisitions, when corn was 
apt to be pre-emergent .or harvested, were not used. There is, however, no way 
of being absolutely sure that these conditions have been avoided. Successful 
generation of crop profile confirmed the imagery growth stage assessment. 
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Acquisitions marred by poor data quality, atmospheric effects, unusual crop 
conditions (such as hail damage or fields cut for silage) or too early or too 
late crop growth stage can be used. The training field will be chosen free 
of these conditions and hence will be unrepresentative of some corn in the 
segment. The corn affected by the conditions will tend to be misclassified; 
this is programmatically correct. 

4.2 ANALYST EVALUATION OF OUTPUT CLASSIFICATION RESULTS 

Input to CLASFY are the line-pixel coordinates of the training field selected 
as a basis for classification, an acquisition set of sufficient temporal dis- 
tribution to characterize the crop profile in the four Landsat channels, and 
initial values for the constants in the modeling function. Outputs are a 
summary sheet of input values, calculated parameters and pixel classification, 
and a full scene classification file, which is translated to a lineprinter 
map of the classification and to a misclassification map of the classification 
file and the ground- truth file. 

• Figure 6 is a copy of the lineprinter summary sheet for sample segment 882. 

• Figure 7 is a copy of the first quadrant of the lineprinter classification 

map. 

a Figure 8 is a copy of the misclassification map of AA pure pixels only. 

• Figure 9 is a copy of the misclassification map of all pixels. 

The classification was evaluated using several criteria. The classification 
map should be "clean." Field patterns should be evident and fields well 
filled out with a minimum of blank spaces (pixels rejected as corn) on the 
interiors. Blank areas should also be clear and be reasonably free of scat- 
tered pixels classified as corn. 

Classification' comparison with the ground truth should have a statistical 
agreement greater than 70 percent and areas of disagreement should be exam- 
ined. In evaluating areas of disagreement, these points are important. 

1. The training field is selected with the constraints that the field should 
not be haze or cloud covered on any acquisition, should not be harvested 
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Figure 9.— Mi sclassifi cation map, all pixels. 
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on the final date, or should not represent an abnormal crop condition 
(such as hail-damaged or cut for silage). Consequently, this method may 
not classify fields that are harvested on the final date or that are ab- 
normal in development into the crop-of-interest category. This is "pro- 
grammatically correct." 

2. Field borders — mixed pixels in general - will tend to be rejected from 
the crop-of-interest category. Again, this is "programmatically correct"; 
the effect is more pronounced in a two-class classifier, such as this, 
than in a maximum likelihood classifier. In this study, by far the 
greatest area of disagreement between the classifications and the ground- 
truth inventory maps was border pixels. Some indication of the number of 
border pixels in an individual segment can be derived from the number of 
AA pure pixels compared to the total number in a Landsat segment; this is, 
however, a subset of the number of pixels which must be defined as mixed 
when multitemporal classification of any type is used. 

Rework techniques were limited to choice of a different set of acquisitions, 
use of four-acquisition classification instead of five, or selection of a 
different training field as the basis of classification. With a satisfactory 
set of acquisitions, training field selection was the most effective rework 
tool. A few classifications were done with relaxed training field criteria 
(this is noted on the tables). The effect of a relaxed chi-square criterion 
seemed to be less distinct field patterns and a more speckled classification 
map. Relaxing the quality-of-the-curve requirements seemed to produce an 
overclassification or an underclassification. 
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5. RESULTS 


The values presented in table 1 must be considered initial results and should 
be viewed as a basis for program development. Evaluation of Badhwar classifi- 
cation as a potential quality assurance check in AA led to the definition of 
a large sample of corn/soybean segments to which this method could be applied. 
As noted in previous sections, within the tentative guidelines, considerable 
experimentation, testing, and exploration of limits and guideline validity 
were conducted. No attempt was made to optimize the procedure based on the 
results. 

Twenty-six of the 29 segments defined for the Multicrop Exploratory Experi- 
ment Test were classified. Results could not be generated for sample seg- 
ments 107 (Boone County, Illinois), 807 (Henry County, Indiana), and 893 
(Webster County, Iowa). 

For each segment, the following items are listed in table 1: 

a'. Sample segment number and location are noted. An asterisk (*) indicates 
that the training fields were defined by the analyst. 

b. Acquisitions available are listed. Consecutive-day acquisitions are 
omitted unless used for classification, and acquisitions outside the grow- 
ing season for corn (Julian days 130 to 300) are omitted. 

c. Acquisitions used for the classification results presented in table l.are 
noted. Landsat-3 acquisitions are denoted by (3). 

d. The coordinates of the training field used for classification and the 
number of pixels in the field are listed. 

e. The percentage of the scene correctly classified (PSCC) based on AA pure 
pixels only is given. 

f. The number of pixels used to compute the percentage above is listed. For 
PSCC based on pure pixels, the base is the number of AA pure pixels in the 
scene minus the number edited by SCREEN minus the number for which there 
is no ground truth. 
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g. A confusion matrix of the classification in scene percentages is given: 

(1) ground-truth corn classified as corn* (2} ground-truth corn classi- 
fied as non-corn, (3) ground-truth non-corn classified as corn, and 

(4) ground-truth non-corn classified as non-corn. These percentages are 
based on calculations by the MISMAP program. 

h. The PSCC based on all pixels, mixed and AA pure, is given. 

i. The number of pixels used to compute the "all pixel" percentages is listed. 
For PSCC "all pixels," the base is 22 932 pixels minus the number edited 

by SCREEN minus the number for which there is no ground truth. 

j. A confusion matrix of the classification of all pixels in scene percentages 
is given. 

k. Additional comments are given in the last column of the table. 
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TABLE 1 CLASSIFICATION RESULTS 


Segment 

and 

location 

Acquisitions 

available 

Acquisition 

used 

Training 

field 

coordinates 

PSCC 

Base 

Confusion matrices, % 




/C - N\ 








(a) 

(a) 

Iffil 

U * n] 

123 

152, 161, 197, 

152(3). 161, 197, 

(3,58), (3,64), 

86.8 

12 897 

24.8 

6.8 

Hamilton, 

233, 269 

233, 269 

(8,66), (8,60), 



6.3 

62.0 

Ind, 



33 pixels 

'82.7 

TsTT? 

24.3 

9.4 







7.7 

58.4 

127* 

152, 161, 197, 

161, 207(3), 

(66,123), (66,133), 

85.9 

16 817 

VHEKi 


Montgomery , 

207, 216, 243, 

216, 243(3), 

(73,136), (73,126), 



iIHEiIwm 


Ind. 

252, 269 

252 

56 pixels 

82.1" 

19 860 

VnTBB 

10.4 







7.2 

41.7 

133* 

152, 197, 233, 

152(3), 197, 

(3,108), (3,115), 

79.3 

16 440 


14.8 

Whitley, 

251, 260, 269 

233, 251 , 269 

(8,117), (8,110), 




64.7 

Ind. 



38 pixels 

75.5 

20 684 


17 7 







6.7 

60.6 

135 

130, 166, 229, 

166, 229, 247(3), 

(25,146), (25,155), 

79.0 

17 599 

20.4 

18.0 

Chickasaw, 

247, 265, 274, 

265(3), 27 4 

(33,156), (33,147), 





iowa 

283, 292 


36 pixels 

75.0 

21 969 

18.7 

nwi 







4.1 

56.3 

141* 

130, 167, 212, 

167, 212(3), 

(14,161), (13,169), 


17 138 

14.8 

7.6 

Madison, 

220, 256, 265, 

220, 265(3), 

(17,172), (18,163), 

■H 


8.9 

68.8 

Iowa 

274, 292 

274 274 

39 pixels 

E M 

19 722 

15.7 

10.8 





m 


9.7 

64.1 

144 

130, 165, 183. 

165, 219, 238, 

(27,166), (27,172), 

90 

15 270 

9.5 

8.3 

Wapello, 

219, 238, 246, 

246(3), 264(3) 

(40,178), (40,171), 



1.3 

80.5 

Iowa 

264, 274, 292, 
300 

39 pixels 

86.4 

17 864 

9.5 

11.6 






1.9 

76.9 

202 

167, 212, 221, 

167, 212(3), 

(17,9), (17,16), 

83.5 



9.7 

Atchison, 

266, 275, 284, 

221, 266(3) 

(20,18), (20,9), 


mm 


72.9 

Mo. 

293 


33 pixels 

80.6 

SSI 

nn 

Ti~4 






■n 

7.0 

69.7 

205 

137, 155, 164, 

155(3), 218, 

(28,182), (28,191), 

84.7 

16 545 


11.4 

Clark* Mo. 

209, 218, 219, 

219, 246(3), 

(33,192), (33,181), 




78.8 


246, 272, ?82, 
290 

272 

37 pixels 

82.5 

18 736 


13.3 

76.0 

216* 

130, 184, 202, 

184, 220, 

(31,184), (31,187). 

87.7 

19 582 

2.6 

KB 

Mercer, Mo. 

220, 238, 247, 

247(3), 274 

(37,187), (37,184), 



9.4 



265, 274, 292 


22 pixels 

86 .6 

20 341 

3.2 








9.7 

83.4 

800 

130, 164, 218, 

164, 219, 

(68,14), (68,34), 

67.2 

16 899 

26.7 

29.7 

Clinton, 

219, 246, 272, 

246(3), 281(3) 

(70,34), (70,15), 



3.0 

40.5 

Iowa 

281, 290, 300 


38 pixels 

65.3 

20 616 

" 24.4 " 

30.8 







3.9 

40.9 

809* 

164, 209, 218, 

164, 218, 

(20,13), (19,23), 

76.8 

14 063 

46.3 

16.4 

Ogle, 111. 

244, 262, 271, 

244(3), 

(27,26), (28,22), 




30.5 


281, 289 

262(3), 271 

26 pixels 

73.8 

17 130 


18.5 







mm 

31.6 

832 

151, 160, 178, 

151(3), 160, 

(91,143), (91,153), 

75.5 

17 335 

12.2 

6.6 

Adams, Ind. 

232, 268, 304 

232, 268 

(94,153), (94,143), 



17.8 

63.3 




39 pixels 

72.5 

20 662 

... 13 s 

8.9 







18.6 

58.6 

837* 

180, 198, 207, 

180, 198, 216, 

(98,35), (96,47), 

80.0 

17 204 

37.5 

6.8 

Benton, 

216, 225, 234, 

234, 252 

(101,48), (104,37), 



13.2 

42.5 

Ind. 

243, 252, 270 


39 pixels 

76.5 

20 409 

35.7 

8.9 







14.5 

40.8 

842* 

160, 176, 232, 

178, 232, 250, 

(40,156), (40,162), 

77.2 

14 436 

29.7 

12.3 

Henry, Ind. 

250, 268 

268 

(45,164), (47,160), 



10.3 

47-5 



32 pixels 

ft** 

17 589 

28.4 

14.9 







12.1 

44.6 


Comments 


Acquisition coverage 
margi nal . Severe 
screening problems. 


Problems with data [ 

quality; additional \ 

early acquisition 
(198) on order. . 


Small Indiana fields; 
misregistration. 


SCREEN removed agri- 
cultural fields; so 
classification is done 
without any screening. 
Marginal acquisition 
coverage and data quality. 
Training field exceeded 
criteria in channel 3 
chi square. 



Screening removed on 
acquisition 165. Poor 
data quality; marginal 
acquisition coverage. 

Acquisitions inadequate. 
Clouds, haze affect 
days 212, 221, and 275. 


Marginal acquisition 
coverage (164 is one-third 
haze, clouds and cannot 
be used). 

Low corn segment, small 
fields, excessive screening 
on days 202 and 265. Chi 
square = 0.0 in channel 1; 
missed some ground-truth 
corn fields. 

Corn crop was hail -damaged 
before acquisition 246; 
considerable harvest evi- 
dent on day 281 . 

Segment more than 5031 corn. 
Training field exceeded chi- 
square criteria in channel 1. 


Inadequate acquisition 
coverage; day 178 is 
cloudy. 


Training field exceeds 
chi-square criteria in 
channels 3 and 4; good 
acquisition coverage. 

Training field exceeds 
chi-square criteria in 
channel 1, more screening 
than apparent data 
quality merits. 


a The PSCC, the base, and the confusion matrices for "pure pixels" appear above the dashed line; for "all pixels," below the dashed line. 
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TABLE 1 Concl uded 


Segment 
! and 

Acquisitions 

available 

Acquisition 

used 

Training 

field 

PSCC 

Base 

Confusion matrices, i 

Comments 

location 

coordinates 



(r -*■ c\ 

/c n\ 






(a) 

(a) 

(5*3 

(5*5) 


843* 

152, 160, 178, 

178, 197, 233, 

(13,121), (13,123), 

80.5 

17 605 

17.6 



Henry, Ind. 

197, 233, 25], 

251 

(22,127), (22,125), 



6.3 




'268 


22 pixels 

77.0 

21 327 

17.4 

15.9 








7.2 

59.6 


852* 

151, 160, 178, 

178, 232, 250, 

(49,36), (49,44), 

82.8 

16 974 

9.1 

15.4 

Training field exceeds 

Randolph, 

232, 250, 268 

268 

(51,46), (53,37), 



1.9 

73.7 

chi-square criteria in 

Ind 



32 pixels 

72.6 

19 492 

1.8 

1.8 

25.6 

70.8 

channels 1 and 2, some 
misregistration, wide range 
of planting dates and much 
of the corn failed to keep 

, 








signature sequence through 
day 178. 

853* 

151, 160, 178, 

160, 178, 232, 

(35,68), (34,72), 

75.6 

15 135 

13.2 

20.7 

Training field exceeds 

Randolph, 

232, 250, 268 

250, 268 

(41,77), (43,71), 



3.4 


criteria in channels 1,2, 

Ind. 



32 pixels 

72.3 

18 552 

13.1 

4.4 

gljtST 

and 4; severe problems with 
chi-square criteria on 







this segment 

860* 

151, 160, 178, 

160, 197, 232, 

(91,61), (91,66), 

82.8 

16 309 

16.7 

11.7 

Misregistration. 

Wells, Ind. 

197, 232, 251, 

251 , 268 

(95,67), (95,62), 



5.5 

66.1 


268 


27 pixels 

78.8 

19 928 

16.8 

14.6 








6.7 

62.0 


864* 

150, 159, 186, 

159(3), 186 

(65,116), (65,120), 

77.5 

15 430 

■fK 



Crawford, 

222, 231, 249, 

222, 231(3), 

(69,120), (69,116), 



iHrXsnl 



Iowa 

258, 267, 294 

267(3) 

24 pixels 

7378 

19 308 

KnrWfil 

22.8 








3.3 

51.4 


865* 

150, 159, 168, 

168, 186, 231(3), 

(65,31), (65,35), 

82.3 

17 562 

22.3 

10.4 


Crawford, 

186, 231, 249, 

249(3), 267(3), 

(72,38), (73,35), 



6.4 

61.0 


Iowa 

267, 294 


32 pixels 

79.0 

21 877 

22.0 

14 0 








7.9 

57 0 


877* 

150, 186, 222, 

186, 222, 231(3), 

(30,61), (28,67), 

72.3 

13 896 

26.1 


Segment has large areas 

Ida, Iowa 

231 , 267 

267(3) 

(33,69), (36,63), 



4.8 


for which ground truth 




38 pixels 

69(3 

17 818 

23.2 

25.1 

was not available. 







5.7 

46.1 


878 

131, 186, 221, 

186, 221, 266(3), 

(9,140), (9,148), 

71.6 

18 360 

18.5 

25.5 

Segment processed as an 

Kossuth, 

266, 293 

293 

(19,151), (19,144), 



3.0 

53.1 

experiment; day 293 too 

I owe 



40 pixels 

7974"" 

21 923 


26.7 

late to be used. 







4.0 

51.4 


880 

150, 186, 204, 

186, 222, 231(3), 

(46,123), (46,133), 

77.4 

mm 

23.2 

21 .8 


Monona, 

222, 231, 249, 

257(3) 

(54,136), (54,127), 


mm 

0.9 

54 2 


l0Wd 

267, 294 

35 pixels 

74.3 


20.9 

24.4 







mm 

1.4 

53.4 


881* 

159, 186, 213, 

159(3), 386, 

(13,81), (12,90), 

83.6 

17 916 

29.1 

14.9 


Monona, 

222, 231, 249, 

222, 231(3), 

(16,91), (17,82), 



1.3 

54.5 


Iowa 

267 

267(3) 

39 pixels 

79.1 

21 740 


18.5 









52.5 


888* 

150, 159, 186, 

159(3), 186, 222, 

(77,77), (76,84), 

87.3 

IS 967 

34.1 

10.7 


Palo Alto, 

213, 222, 231, 

231(3), 267(3) 

(81,87), (82,80), 




53.2 


Iowa 

258, 267, 293 

36 pixels 

83.3 

19 057 


13.6 








1H 

51 .4 


891 

15?, 168, 186, 

168, 186, 204, 

(90,5), (90,13), 

73.4 

14 828 

25.7 

flBVjjfe 

Training fields exceeded 

Shelby, 

204, 249, 258, 

258, 267(3) 

(97,13), (97,5), 



4.5 


chi-square criteria in 

Iowa 

267 


37 pixels 

70.0 

19 010 

23.9 

5.9 

US 

channels 3 and 4; data 
quality marginal . 


a The PSCC, the base, and the confusion matrices for “pure pixels" appear above the dashed line; for "all pixels," below the dashed line. 
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6. CONCLUSIONS AND RECOMMENDATIONS 


Badhwar classification provides an approach to classification different 
from that of the classification methods currently in use. Basic data require- 
ments are similar to those of current methods, however; and Badhwar classifi- 
cation is applicable to a very large percentage of LACIE sample segments. 

This paper documents initial results as a base for program development. Accu- 
racy averaged over the 26 segments is 80.4 percent for AA pure pixels only 
and 76.6 percent for all pixels. 

Further research using Badhwar classification should be conducted on a subset 
of' the site data set used for this study. The following set of sample seg- 
ments is recommended: 


127 

852 

133 

853 

135 

860 

141 

864 

205 

865 

216 

877 

809 

880 

837 

881 

842 

882 

843 



The available acquisitions for these particular sites seem to be well suited 
for meeting the criteria of data quality, acquisition distribution, and 
coverage of the corn growth cycle. 

Further evaluation of this classification technique is recommended. 
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